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Abstract 

In this note we further investigate the procedure for computing tree-level amplitudes in 
Yang-Mills theory from connected instantons in the B-model on P 3 ! 4 , emphasizing that the 
problem of calculating Feynman diagrams is recast into the problem of finding solutions to 
a certain set of algebraic equations. We show that the B-model correctly reproduces all 6- 
particle amplitudes, including non-MHV amplitudes with three negative and three positive 
helicity gluons. As a further check, we also show that n-particle amplitudes obtained from 
the B-model obey a number of properties required of gauge theory, such as parity symmetry 
(which relates an integral over degree d curves to one over degree n — d — 2 curves) and 
the soft and collinear gluon poles. 
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1. Introduction 

In Witten proposed a remarkable connection between scattering amplitudes in 
Yang-Mills (YM) theory and a certain topological string theory, the B-model on P 3 ' 4 (recent 
related work includes [0-0). This conjecture leads to the following formula (equivalent to 
one first written down in and studied further in 0,g|) for the n-particle amplitude, 
written in a manifestly Af = 4 supersymmetric notation: 

n— 3 „ n d / n \ / n \ 

A n = z(2rr) 4 ^ M 2 £ / dM n , d ]J 5 2 (Xf - ^If) ]J S 2 ^ fc^A? U 4 ^ • 

d=l ^ i=l fc=0 \i=l / \i=l / 

(1.1) 

The details of this formula will be clarified in section 2, but we have written it down 
here in order to stress its simplicity and importance. We believe that (|1 . 1| ) encapsulates 
the complete n-particle tree-level S'-matrix of YM theory (for any gauge group), thereby 
providing an exact solution of classical YM theory in four dimensions. This formula sums 
up a huge number of Feynman diagrams (see for example Fig. 1) into an expression which 
fits on a single line. In this note we provide strong evidence supporting our confidence in 
this formula and explore some of its structure. 

The formula ( |1 . 1|) was derived by considering the contribution to the scattering am- 
plitude from a single connected instanton (= holomorphic curve in P 3 ! 4 ) of degree d in the 
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A 6 = ^xr^xj-i Lrv^rxr + jvwwWuww + 218 diagrams 



Fig. 1: The standard computation of a six-gluon tree amplitude requires 
summing 220 Feynman diagrams (in conventional gauges) . 



topological B-model. (As explained in JIJ , counting the fermionic modes reveals that the de- 
gree d is related to the total helicity hi °f an n particles according to d = \{n— Y hi— 2).) 
In |I] Witten speculated that one might have to consider, in addition to (|1.1|), contribu- 
tions from collections of disconnected instantons of degrees d\ with J^di = d. (See Fig. 2 
for a schematic depiction for the n = 6, d = 3 amplitude.) 

However, it was found in HQ that the formula ( |1 . 1|) correctly reproduces the known 
YM result for the mostly minus MHV (maximally-helicity violating) amplitudes (some- 
times called 'googly') that are related to the mostly plus MHV amplitudes by complex 
conjugation. Even though the googly amplitudes are calculated from an integral over the 
moduli space of instantons of arbitrarily high degree, precise agreement was found with 
gauge theory without the need for additional contributions. 

More recently, a novel method for calculating YM tree amplitudes, also motivated by 
the B-model on P 3 ! 4 , was proposed in a very interesting paper ||. The starting point for 
their proposal involved considering only completely disconnected instantons (i.e., d instan- 
tons of degree 1). Remarkably, it was found that their rule also gives correct gauge theory 
amplitudes. The B-model seems to give two separately correct methods for calculating YM 
tree amplitudes, rather than a set of contributions which need to be summed (see Fig. 2). 

The proposal of |J involves a diagrammatic expansion which bears no apparently obvi- 
ous connection to the formula (|1.1| ), except that they both seem to be correct. It would be 
very interesting to understand directly the relation between these two methods. Moreover, 
if the B-model on P 3 ! 4 gives us two not obviously equivalent formulas for YM amplitudes, 
then it will likely give us an infinite family of formulas (which roughly speaking weight the 
different types of diagrams in Fig. 2 differently). Undoubtedly we have only encountered 
the tip of the iceberg connecting the topological B-model to Yang-Mills amplitudes. 

So far the formula ( |1.1|) had only been checked for MHV and googly amplitudes P,§| . 
In section 3 of this paper we confirm that the formula also gives the correct 6-particle non- 
MHV amplitudes. In section 4 we check that for any n and d, (|1 . 1| ) satisfies a number of 
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Fig. 2: Schematic depiction of how one might have thought to organize the 
calculation of the 6-particle mostly minus MHV (d = 3) amplitude in the B- 
model on P 3 ' 4 . The dark x's mark the insertions of the 6 external particles, 
the dotted line is a twistor space propagator (constructed in ||), and the 
solid lines represent instantons (i.e., holomorphic curves in P 3 ! 4 ) of degree 
d = 1, 2, 3 (schematically encoded in the waviness of the curve). Although one 
might have expected that it would be necessary to sum together contributions 
of all three types, we find that the single diagram of the first type (studied 
here and in HH), and the sum of the 21 diagrams of the third type (studied 
in 0), separately give the correct gauge theory answer. 

properties required of general Yang-Mills amplitudes, such as the soft and collinear gluon 
limits. Of particular importance is parity symmetry, which requires that (|1 . 1] ) should be 
invariant under A <-> A. This non-manifest symmetry of (|1.1| ) is proven explicitly in section 
4.2 below. We conclude with a list of open questions and puzzles. First, however, we turn 
our attention to the details of ( p. . 1|) and highlight a crucial fact about the formula: namely, 
that it is not really an integral at all. 



2. The Main Formula 

In this section we first clarify the ingredients appearing in the formula (|1 . 1| ) and 
then investigate some of its mathematical properties. The quantity A n in ( |1 . 1| ) denotes 
the color-stripped n-particle partial amplitude (see for example [PI), and we employ the 



spinor helicity notation in writing A n as a function of (\f, \f, ^a)) i = 1, ■ ■ - ,n, where 
A and A are commuting real two-component spinors of positive and negative chirality, 
respectivelyS, and tja is the four-component Grassmann coordinate of Af = 4 superspace. 
The P a are two degree d polynomials in a which we parametrize as 

d 



1 For simplicity we work in signature + H , where such spinors are possible. It is straight- 
forward to analytically continue the tree-level YM amplitudes to signature H if desired. 
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in terms of 2d + 2 coefficients (moduli) a£. When needed, we will follow the conventions 
of @U in denoting P} = A, and P? = B t . The measure for integration in ( |1 . 1|) is 

dM n , d = unT7oY\ 11 F? V 2 " 2 

vol(GL(2)) f^&to - £7 i+ i) 

The factor of l/vol(GL(2)) is included because the integrand is invariant under a certain 
GL(2) symmetry and so the integral would otherwise be infinite. Practically, the conse- 
quence of this factor is simply that we can choose to fix four of the variables (any one of 
the a's and any three of the cr's) at the expense of introducing a Jacobian factor of 

J = o(crj - cfj)(cfj - a k )(cf k - o-j). (2.3) 

The choice of which a and which three a's to leave un-integrated is arbitrary and does not 
affect the final result. 

2.1. A key point 

The single most important fact about the integral ( |1 . 1| ) is that it is not really an 
integral. To see this, let us start by showing that (|1 . 1] ) respects momentum conservation. 
Taking a particular linear combination of the quantities set to zero by the delta functions 
gives 

d / n \ n n n 

= E a * = = E ^ = E^ Q , (2.4) 

k=0 \t=l / i=l i=l i=l 

where we used the definition ( |2.1] ) and some more delta functions from ( p. . 1|) . Therefore, 
the delta functions in ( |1.1|) indeed force overall momentum conservation. 

At the practical level, this means we can 'pull out' the overall factor of S 4 (^2pi) at 
the expense of introducing a Jacobian, by using an identity such as 



% "i t 



i=l v 1 l/ fc=0 \i=l 1 



A? BA A.o/AAfA, 1 ^ 



(2.5) 



,i=l / j=3 v 1 l/ k=l \i=l 



(where we used & = Xj/Ai). In writing this identity we have made a particular choice 
of which four delta functions to pull out. There is however no canonical choice, and 
different choices are useful for different calculations (and lead to different Jacobians), so it 
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is convenient to leave momentum conservation slightly scrambled into the delta functions 
in Note that supermomentum conservation 5 8 (Yl ^hVia) pulls out similarly. 

Let us now return to the claim that ( |1 . 1] ) is not really an integral. The measure dJ\A nd 
in fl2.2p has (2d+ 2) + (n) + (n) — (4) = 2n + 2d — 2 integration variables, while the integrand 
in ( |1 . 1|) has 2n + 2d + 2 delta functions. If we 'pull out' the overall momentum conservation 
delta functions, then for any n and d there are precisely as many integration variables as 
delta functions. Therefore the entire integral is supported on a discrete set of points, and 
the formula fll.lj ) is just a recipe to solve the 2n + 2d + 2 polynomial equations 

d 

\° = &X>fa*, a = 1 ' 2 ' i = 1 ' 



,n, 

(2.6) 

= ^&ofAf, a = 1,2, k = 0,...d 



k=0 

n 



i=l 



for the variables (a^,cri,^i), and then to sum a certain Jacobian (obtained in the usual 
way from (|1 . 1| ) ) over the collection of roots. 

One of the most interesting questions about the system (|2.6|) is: what is the number 



of roots N n> d for general n and dl At this point all we know for sure is that 

N n>1 = iV n>n _ 3 = 1, iV 6j2 = 4. (2.7) 

The first two cases are MHV and googly amplitudes previously studied in the literature, 
and iV 6j 2 is the non-MHV 6-particle amplitude discussed in the following section. In section 
4 we prove that N n ^ n -d-2 = N n ^- Certainly it would be very interesting to have a better 
understanding of the mathematics underlying the equations ( |2.6| ). In particular, it would 
be especially interesting to learn how N n> d grows with n and d. 

2.2. A complex puzzle 

A priori, the moduli of the curve and the coordinates o~i on P 1 should all be complex 
variables. In order to evaluate the integral ( |1 . 1| ) it is necessary to specify an integration 

contour in this 2n + 2d — 2 complex dimensional space. In spacetime signature + H it 

makes sense to take A and A to be independent real variables, and it is natural to choose 
the integration contour in which all of the a£ and o~i are real. 

For both the MHV [d = 1) and googly (d = n — 3) cases, the unique root of the 
equations ( |2.6|) indeed has the property that o~i and a% are real. However, for the 6- 
particle amplitude with d = 2, which we discuss in section 3, there is a puzzle. Depending 
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on the choice of A and A, there can be four real roots, two real roots and one complex 
conjugate pair, or two complex conjugate pairs. The YM tree amplitude, which is always 
real (forgetting the % in front of fll.l] )), is reproduced only if all four roots are summed 
over, regardless of whether they are real or complex. 

The lesson from this analysis is that restricting ( |1 . 1| ) to the contour where all a's and 
a's are real does not give the correct gauge theory scattering amplitudes. In fact, we do not 
know how to write any contour which makes the integral formula ( |1 . 1| ) valid for arbitrary 
choices of A and A. This amplifies the comment we made at the beginning of the previous 
subsection: the formula ([T[l]) is n °t really an integral. To overcome this problem we avoid 
thinking about ( |1 . 1| ) as an honest integral, but instead view it as a recipe for finding the 



solution (which in general can be complex) of (|2.6|) and then summing a Jacobian over the 
set of roots. 

2.3. A diagrammatic expansion? 

From this new standpoint, let us ask ourselves whether the formula might have an- 
other, more natural interpretation. The fact that the computation of a scattering am- 
plitude from the formula (|1.1|) reduces to summing a certain quantity over a finite set of 



points is reminiscent of some sort of diagrammatic expansion, where, for example, ( |2.7| ) 
suggests that there is a single diagram for mostly plus and mostly minus MHV amplitudes, 
while four diagrams contribute to the 6-point non-MHV amplitudes. 

It is tempting to wonder whether there is any connection between such 'diagrams' 
and the new diagrammatic expansion for YM scattering amplitudes which was recently 
proposed in . According to their proposal, A n ^ is associated with the collection of trees 
with n cyclically labeled external legs and d vertices, such that each vertex has at least 3 
legs. For general n and d there are ^ ( jZi) { n ~d-i 2 ) such graphs!, which in all cases except 
the trivial case d = 1 is larger than (|2.7| ). (We have written the number of diagrams in 
M = 4 superspace. For particular choices of helicities of the external particles there are 
frequently fewer diagrams.) 

However, the diagrams of || have an additional symmetry in the form of an arbitrary 
spinor rj a which drops out only after summing together all of the graphs. The number of 



The counting of these graphs is equivalent to a combinatorial problem which appeared in 



Plutarch's biographical notes on Hipparchus [12]. We are grateful to C. Herzog for many fun and 



enlightening discussions regarding the combinatorics. 
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diagrams is not gauge invariant, and special choices of r\ can set whole classes of diagrams 
to zero. In contrast, our 'diagrams' have no residual manifest symmetry — the GL(2) 
cancels out diagram by diagram (root by root) and does not change their number. Maybe 
there is some choice of r\ for which the diagrams of || reduce, in number and in value, to 
the contributions obtained from the roots of our formula (|1 . 1[) ? 

We believe it is more likely that the topological B-model has some huge symmetry 
group which relates the formula ( |1 . 1| ) , with its associated 'diagrams', to the diagrammatic 
expansion of |J. Their parameter r\ is a small residue of that huge symmetry. 



3. The 6-Particle Non-MHV Amplitudes 

In the previous section we introduced the formula and discussed its basic proper- 
ties. But what is the connection between ( |1 . 1| ) and the n-particle scattering amplitude in 



gauge theory? In M it was shown that a prescription equivalent to the d = 1 case of (1.1 



reproduces the mostly plus MHV amplitudes in YM theory. In it was shown that the 
formula also works for mostly minus MHV amplitudes (sometimes called googly or MHV). 
These have d = n — 3 and are related (in Minkowski signature) to MHV amplitudes by 
complex conjugation. 

Although the latter check involved an apparently nontrivial integral over the moduli 
space of curves of arbitrary degree in P 3 ' 4 , the question of whether ( |1 . 1| ) is correct for 
genuinely non-MHV amplitudes was left open. The simplest amplitudes which are neither 
MHV nor googly are those with n = 6 particles and d = 2. Since we work in a manifestly 
M = 4 formalism, our results apply simultaneously to all possible helicity orderings (when 

all six particles are gluons, there are three independent helicity orderings: + + H , 

+ H 1 and + - + - + -). 

In this paper we report that the formula (|1 . 1| ) , in the case n = 6 and d = 2, precisely 
matches the 6-gluon scattering amplitudes first computed by Mangano, Parke and Xu 
||13|| . We originally obtained this result numerically, by (1) choosing at random a collection 
of (Aj, Aj) (subject to overall momentum conservation (|2.4|) ), (2) numerically solving the 
polynomial equations ( |2.6|) , which were always observed to have four roots, and then (3) 
summing the Jacobian obtained from (|1 . 1|) over the four roots. The whole calculation takes 
only a few seconds on a fast computer and can be repeated as often as desired for different 



(Aj, Aj). The result was always found to agree spectacularly with the formula given in [|T3|] . 
Note that all three independent helicity configurations can be checked at the same time 
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since the choice of helicities only affects the fermion determinant and does not change the 
value of the roots. 

The only puzzle we encountered is that occasionally, for some (Aj, Aj), the roots are 
complex, as we discussed in subsection 2.2. Precise agreement with gauge theory was 
nevertheless always found by doing the most naive thing possible and summing over all 
four roots, whether real or complex. 

Unfortunately, it seems rather difficult to construct an analytic proof that the formula 
( |1 . 1|) is correct for the case n = 6, d = 2. Let us now outline the best line of attack that 
we know of at the moment. We will not give precise formulas for each intermediate step 
because they are extremely lengthy and moreover because we are hopeful that a more 
clever way of analyzing the equations will become available. We believe that only after the 
mathematical structure of the equations (|2.6| ) is better understood (for arbitrary n and d) 
will it be clear how best to organize this calculation analytically. 

3.1. Constructing a Groebner basis: a sketch 

The most interesting result of the numerical analysis is that the number of roots is 
-^6,2 = 4, which does not appear obvious from ( |2.6|) . Recall that we can fix one of the 
a's and three of the cr's (say o~i, 02 and 03) using the GL(2) symmetry. The remaining 
In + 2d — 2 = 14 'integration variables' are fixed by solving ( |2.6| ). In fact, it turns out 
to be possible to express all of the a's, £'s and two of the remaining three a's as rational 
functions of the final a (say Moreover, one can extract from (|2.6|) a single equation 
which is quartic in ct 6 and does not depend on any of the other 'integration variables.' 
The coefficients of this quartic polynomial are themselves polynomials in a%, o~2 and 03 
and the covariant kinematic quantities [ij] and (ij). The four roots of this master quartic 
equation determine the solutions for all 14 variables. 

Here is a schematic description of how to derive this quartic equation. Choose some 
subset S of the equations and use them to solve for the a's and £'s in terms of the a's. 
Plugging the solution into the remaining equations gives polynomial equations just on the 
a's. This process can be repeated many times by starting with different sets S of equations, 
leading to a large number of polynomial equations on o~4, a$ and ctq. The game then is to 
find the common roots of these polynomial equations. In mathematical language, we need 
to construct a Groebner basis for the ideal generated by these polynomials. Let us now be 
a little more specific. 
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Start with the equations on the top line of (|2.6| ) . By eliminating £j between the a = 1 
and a = 2 versions of this equation, one arrives at the six equations 

2 2 
k=0 k=0 

which are conveniently expressed in matrix notation as 

f\\ \\a x \\a\ A? \\a x A?a?N 

\1 \1_ \l r ,2 \2 \2_ \2^2 

A 2 A 2 a 2 A 2 a 2 ^2 A 2 a 2 A 2 a 2 

A 3 A3O3 Agaf A 2 A 2 a 2 A 2 a| 

A 4 A 4 a 4 A 4 a 4 A 4 A 4 cr 4 A 4 a 4 

A5 Xla 5 A 5 cr 5 A 5 A|a 5 A 5 cr| 

\Ag Aga 6 Aga| A| A|cr 6 

A nontrivial solution exists if and only if the determinant of this matrix is zero: 

= X= ^ ^ijklmnV(iJ,kJ,m,n){il){jm)(kn). (3.3) 

i,j,k,l,m,n 

Here V is the cyclic product of cr's (not the Vandermonde matrix), 

V(i, j, k, l,m,n) = (cr* - <7j)((7j - cr fc )(cr fc - cri)(cri - a m )(cr m ~ &n)(cTn - <?i)- (3.4) 

Another way to think about this equation is as follows. Since one of the a's is fixed by the 
GL(2) symmetry, we really only are allowed to solve for five of the a's. If we choose any 
five of the equations (|0]) to solve for the five a's and then plug the solution into the sixth 
equation, we find the condition that ( |3.4| ) should vanish. 

Next we turn our attention to the equations on the second line of ( |2.6[ ) . These are six 
(a = 1, 2, k = 0, 1, 2) homogeneous linear equations on the six variables £j. When cast in 
matrix form, the relevant matrix is precisely the transpose of (|3.2| ), but with A «-> A. A 
nontrivial solution exists if and only if the corresponding determinant vanishes: 

= X= ^ €ijkimnV{i,j,k,l,m,n)[i j)[kl)[mn). (3.5) 

i,j,k,l,m,n 

So far we have obtained (subject to X = 0) a unique solution for all of the moduli a£, 
and (subject to X = 0) a unique solution for all of the The final step is to require that 
these solutions are compatible, in that they obey the top line of (|2.6|) . There are a huge 
number of such compatibility conditions that one can form, depending on which five of the 
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six equations (|3.1] ) one uses to solve for the moduli and which five of the six equations from 
the second line of ( |2.6| ) that one uses to solve for the £j. These equations are polynomials 
in 04, cr 5 and 0% whose coefficients depend on A, A and the fixed values of 01, 02 and 03. 

However, these equations (as well as the X = = X equations) all have spurious 
roots at (T4 = (T5 = o"6 = cr, for i = 1,2, 3. To eliminate these roots one constructs a linear 
combination of these equations (with coefficients involving powers of 04 and 05), with the 
coefficients chosen so that the result factors into a single quartic polynomial q , (cr 6 ) without 
the spurious roots times a high-degree polynomial with only spurious roots. 

In the previous few paragraphs we have explained in words the process of constructing 
a Groebner basis for the ideal generated by the polynomials ( [2.6| ). Once the roots are 
found, it remains to evaluate the Jacobian. At the end of the day, the amplitude can be 
written schematically as a rational function in erg, summed over ag satisfying some quartic 
polynomial: 

Ae,,= £ 44 (3.6) 

{a 6 :q(a 6 )=0} 

Abel's theorem guarantees that the result of this sum is a rational function of the coef- 
ficients of the polynomials p, q and r, and it is easy to check numerically that the result 
precisely matches the gauge theory amplitude of Mangano, Parke and Xu fl3]. More gen- 



erally, Abel's theorem guarantees that for any n and d, (|1 . 1|) turns into a rational function 
of the covariant quantities (ij) and [ij] once all of the roots of (|2.6|) are summed over. 

3.2. Analysis for special A 

Although the n = 6, d = 2 amplitude is complicated in general, instructive analytic 
expressions can be obtained by considering special cases. For example, let us here consider 
the case Xf = A| = and [1 5] = 0. For this degenerate case, numerical investigation 
reveals that there are only three roots (one is a double root — the statement that N^2 = 4 
is always true when one counts multiplicities) . Let us demonstrate analytically how to find 
these three roots. 

We fix the GL(2) symmetry by setting a\ = 1 and a\ = {0,1,-1} for i = 1,2,3. 
Also, without loss of generality we can rescale the A's to set A| = 1. From the AiXi = Bi 
equations for i = 2, 3, 4, 5, 6 we can solve for the moduli a\, a\, af, o| an d 04 in terms of 
05, o"6. The first solution is a 4 = 0, and the other one is 

^ _ A 3 A 5 A 6 cr 53 cr 6 3056 + A2(2A3A 6 crfcr 6 2cr 6 3-A 5 cr52(A 6 a5 6 q- 6 3 + 2X 3 a 53 al)) 

A3A5A6<753a" 6 30"56 + A2(A5A60"520"62C56 + 2As(— (7^X^52^53 + A6Cr50"62C63) 
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The 04 = root gives a unique solution for <r 5 , when we plug the expressions for 
a\,a\,a\ and a| into the equations following from the second line of (|2.6|) . The nonzero 
04 root gives a simple solution for cr 6 : 

= [65](A 3 -A 2 )A 6 

6 2[45]A 2 A3 + [65](2A 2 A 3 -A 2 A 6 -A3A 6 )' 1 ' ] 

and a quadratic equation on 05. In other words, the analog of the fourth order polynomial 
described in the previous subsection factorizes into a quadratic one and the square of a 
linear one. Solving the equations and plugging them into the Jacobian gives a result which 
agrees numerically with the known gauge theory result. 



4. Checks on n-Particle Amplitudes 

To summarize, we now know that the formula ( |1 . 1|) correctly reproduces all MHV 



and MHV amplitudes, as well as all 6-particle amplitudes. The nontriviality of these 
checks makes it implausible that some complication arises for further amplitudes which 
might render ( |1 . 1|) invalid. Nevertheless, it would certainly be satisfying to prove that 
the formula ( |1 . 1| ) is correct, perhaps by showing that it satisfies the recursion relation of 
TM. Since we do not have a complete proof yet, we will content ourselves with tabulating 



several consistency checks that ( p..l| ) is indeed the tree-level S'-matrix of YM theory for 
arbitrary n and d. 

4-1. Some properties of gauge theory scattering amplitudes 

Color-ordered scattering amplitudes in YM theory satisfy a number of important 
properties, including: 

(i) Cyclicity: 

A(2,3,...,n,l) = A(l,2,...,n). (4.1) 

(ii) Reflection: 

4(n,n-l,...,l) = (-l)M(l,2,...,n). (4.2) 

(iii) Conjugation: Parity symmetry implies that the amplitude is invariant under inter- 
changing each helicity + «-> — and simultaneously interchanging A <-> A,!. The M = 4 
supersymmetric version of this statement is 



A(X i7 Xi,r]iA) = J d 4n tp exp 



n 

I 



A(Xi, Xi, tpf). (4.3) 



3 This transformation makes sense with our choice of signature (see footnote 1). In Minkowski 
signature the left- and right-hand sides would be related by complex conjugation. 
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(iv) Dual Ward (or Sub-Cyclic) Identity: 

A(l,2,3,...,n) = 0, (4.4) 

C(l,...,n-1) 

where n is held fixed in the last position and C(l, . . . ,n — 1) denotes the set of cyclic 
permutations of {1, . . . , n — 1}. This identity expresses decoupling of the U(l) degree of 



freedom [15 



(v) In [|l6j it was proven that YM amplitudes satisfy the following generalization of (iv): 

^2 A(i 1 ,...,i m J 1 ,...J k ,n + l) = 0, l<m<n-l, m + k = n, (4.5) 

Perm(i,j) 

where the sum is taken over permutations of the set (ii, . . . ,i m , ji, . . . , jk) which preserve 
the order of the (zi, . . . , i m ) and (J±, . . . ,jk) separately. 

(vi) Soft-Gluon Limit: In the limit p\ — > 0, the amplitude behaves as 

(n 2) 

A(l+, 2, . . . , n) — ( ^ 1)(1 ; 2) A(2, . . . , n). (4.6) 

Of course a conjugated version of this equation should also hold in the case when particle 
1 has negative helicity. We do not consider that case directly in this paper, since it follows 
as a result of (iii) above. 

(vii) Collinear Limit: In the limit pi — > zp and p2 — > (1 — z)p for z G (0, 1) and some p 
with p 2 = 0, the amplitude behaves as 

A(l+ 2+, 3, . . . , n) — ^=L=-i-A(p+, 3, . . . , n). (4.7) 

Again it follows from (iii) that there is an obvious conjugate to this relation for the case 
when particles 1 and 2 both have negative helicity. The final case, when particles 1 and 2 
have opposite helicity, is 

z 2 1 (1 — z) 2 1 

A(l + , 2", 3, . . . , n) — ► - — ;A(p+, 3, . . . , n) + \ ' -, — rA(p~,3, . . . , n). 
1 J ^(1 - z) [1 2] ^ J V /I(T37)(12> ^ J 

(4.8) 

(viii) Multi-particle Poles: Color-ordered amplitudes can only have poles in channels corre- 
sponding to a sum of cyclically adjacent momenta going on-shell [[TTJ. That is, if we denote 
Pl,m = Pi + P2 + ■ ■ ■ + Pmi then the amplitude factors in the p\ m — >■ limit according to 

An(l,...,n) — > ^2 ^m+i(l, • • -,m,p x )^— A n _ m+l {m + 1, . . . ,n,p~ x ). (4.9) 
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Properties (i), (ii), (iv) and (v) are manifest in (LI) due to the way the Oi enter in 



( |2.2|) . Indeed they follow so trivially from (|2.2|) that the reader may well wonder why we 
have bothered to mention them. We have done so only because not all of these properties 
are immediately obvious from the Feynman diagram expansion of gauge theory amplitudes. 
(These properties are also not all manifest in the diagrammatic prescription of 

Of the remaining properties, (iii), (vi) and (vii) will be proven in the following sub- 
sections. The final property (viii) regarding multi-particle poles will not be addressed 
here. Indeed, note that a proof that (|1 . 1| ) satisfies (viii) would essentially be a proof that 
( |L1| ) is correct — since a tree-level YM amplitude is uniquely fixed by its poles (and their 
residues) . 



4-2. Parity symmetry 



The parity symmetry ( fL3| ) is obvious in gauge theory but not manifest in the formula 
( |LT|) I. On the individual component amplitudes A n ^, ( [13|) says that 



A n!d (X,X,r]) = I d 4n ip exp 



A-n^n — d— 2 



(4.10) 



thereby relating an integral over the moduli space of degree d curves to an integral over 
the moduli space of degree n — d — 2 curves. 

The proof of ( [I.10|) is fairly straightforward. We start by looking for a way to relate 
the bosonic part of the amplitudes, 



/n d / n 

i=l k=0 \i=l 

/n n—d—2 / n 

dM n , n -d-2l[6 2 (xt -£?*) n 5 2 [y^la\\t ) . 

i=l 1=0 \i=l 



(4.11) 



Here dAi n -d-2 and P are the obvious generalizations of ( |2T2| ) and ( p.l|) : 



d 2(n-d-2)+2~ d n~ d n^n 
dM n -d-2 = ,^ r ^ -[[ 



1 



n-d-2 



ur<rr-v^ II ~ P t = > ■ (4-12) 

vol(GL(2)) t\ U 



The parity symmetry was also very recently discussed in [|| in the framework of 
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We will show that after integrating out the moduli a, the first set of delta functions in A 
exactly transform into the second set of delta functions in A (and vice versa) when one 
makes the change of variables 

0^ = <Ti, li = — 7- (4.13) 

The Jacobian for this coordinate transformation is unity, but we will pick up a simple 
Jacobian from manipulating the bosonic delta functions. This Jacobian will exactly cancel 
a similar fermionic determinant. 

Let us begin by studying the quantity 

We claim that p m is a polynomial in the aVs of degree m — n + 1. To see this, consider p m 
as an analytic function of z = a n (this can of course be repeated for all of the a's). It looks 
like Pm{z) might have poles at the other a i: but in fact it is easy to see that the residue 
is always zero. So Pm{z) has no poles, and grows at infinity like z m ~ n+1 , so it must be a 
polynomial of degree m — n + 1. In particular, p m vanishes for m < n — 1, and p n -i = 1. 

Now consider the first type of delta function in A, (we focus on one value of a and 
restore covariance later) 

/n 
d d+1 al]^-^), (4-15) 

i=l 

and take linear combinations of the delta functions according to the n x n matrix with 
entries 



a™ 



M mi = —= f -, i = l,...,n, m = 0, ...,n- 1. (4.16) 

£'•11.,,,^' - °i) 



That is, we write 



/n— 1 / n \ 
d d+1 a (detM) J] 5 ^M mi (A, - . (4.17) 

m=0 \i=l / 



The second term in the delta function is now 



S p n h — ^ S ak ° k - = E a kPk+m, (4.18) 
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using the definitions (p.l|), ( [4.14j) and ( |4.16| ). Then recalling that Pk+ m is zero for m < 
n — d — 1, we can split the delta functions into two kinds: 

n—d—2 / n \ r n — 1 / n d 

I = (det M) H 8 IJ2 &T X i / dd+la II ME S^Am - E a ^+' 

m=0 \i=l / J m=n-d-l \i=l k=0 

(4.19) 

The d+1 moduli now appear linearly in the last d+1 delta functions and can be integrated 
out trivially. The Jacobian for this is just 1, because Pk+m is a triangular matrix with 
diagonal entries p n -i = 1- Finally, we conclude that 



/n 
i=i 



i=l 



-d-2 



n mes^ 



(4.20) 



m=0 \i=l 



where V is the Vandermonde determinant of all of the cr's and the term in brackets comes 
from evaluating det(M). 

The next step is to simply apply (|4.20|) in reverse to get 

-i 



n* e^^ 

k=0 \i=l 



i=l 



/I Ij 
d^-^j+igjj^-SPO- ( 4 -21) 
i=i 



Finally we can combine ( f4.20|) and ( 4.21 ) and restore the a and a indices to arrive at 



/n d / n 

dM n4 i[5 2 (xt-^pr)U 62 

i=l k=0 \i=l 



dM 



n,n — d—2 



-I 2 



n 



n-d-2 



(4.22) 



z=o 



,i=l 



We have now related the bosonic integral over degree d curves to the bosonic integral 
over degree n — d — 2 curves, up to a factor which with the help of ([4 . 1 3| ) can be written as 



nl 

i=i ? * 



(4.23) 



In fact, this is precisely the factor which should arise from the fermionic Fourier transform 
in the formula G4.10| ): 



d An il; exp 



i=l 



-d-2 



n 5 * 



£=0 



.1=1 



i=i 



n* 4 e^^ 



fc=0 \i=l 



(4.24) 
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This completes the proof that satisfies the conjugation property ( f4.10|) . 

Incidentally, the above arguments show that given any solution of the equations ( |2.6| ) 
one can construct a solution of the conjugate equations 

n-d-2 

Af = 5 a = 1,2, i = l,...,n, 

1=0 (4.25) 

= J^^Af, a = 1,2, l = 0,...,n-d-2 

i=l 



by taking a and £ to be given by (|4.13|) . It is not necessary to independently specify the 



af since the top equations in ( |4. 25|) determine them uniquely in terms of (<7j,£j). Thus, 
we have shown that 

iV nin _ d _ 2 = N n4 , (4.26) 

and moreover, that the contribution to from any given root is exactly the conjugate 
of the contribution of that root to A nn _d-2- The relation (4.26) is reminiscent of the 



relation between Betti numbers for a manifold of dimension n — 2 as well as of the relation 
between Hodge numbers under mirror symmetry. It would be interesting to find a relation 
between N Uy d and some invariants of P 3 ' 4 (perhaps Gromov-Witten invariants) or of its 
moduli space of holomorphic curves. 

4-3. The soft gluon limit 

For both the soft gluon and collinear limits, comparing the left- and right-hand sides 



of (|4.6| ) and (|4.7| ) reveals that we will have to perform two integrals and eliminate two 
delta functions. Clearly we want to eliminate the appearance of gluon number 1 on the 
right-hand side, so we should eliminate the two delta functions 8 2 (Xf — ^l-PfO (f° r a = 1, 2) 
by performing the integrals over £i and o~\. In general there can be several roots which 
contribute to this integral. However, we are only interested in roots which in the desired 
limit give rise to a pole in the amplitude. We will argue that only one root contributes to 
the coefficient of this pole. 

A prototype for both the soft gluon and collinear limits involves an integral of the 

form 

(ii)^oJ a! -ai JK l> \X\Xj [A x Ail) 
Specifically, we are interested in the poles of this integral. We do not yet need the explicit 
form of A, B or /, and need only to make assumptions which are completely reasonable 
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for the application at hand: B/A is a rational function of a with isolated roots, and the 
function / has no poles in a\. 

The quantity in brackets in Q4.27|) vanishes when o\ = cr^ and hence can be written as 



Bi 
Ax 



Bi 



(ai - a i )F(a 1 



for some F. Changing integration variables from o\ to w = a% — ai gives 



lim 

<li)->o 



w 



9{w) 



(jl) 
A1A1 



— wF(w, ai 



(4.28) 



(4.29) 



In the limit (lz) — > the roots of g(w) are easy to analyze. There is one root (which we 
will call w = wq) for which w is small (of the same order as (1 £)), and there may be other 
roots for which F(w, ai) is small. We assume there is no degeneracy amongst the possible 
roots. Integrating the delta function gives a factor of l/g'(w), which is a number of order 
1 at any of the roots. Therefore, the only pole in the integral Ii comes from the factor of 
1/w evaluated on the root w = wo — > 0. 

The value of wo is given by the implicit equation 

(il) 1 



w 



\{\}F{w ,<Tiy 

with F(wo,ai) being of order unity. The contribution of this root to the integral is 



(4.30) 



/ = 1 (?± 

ow 



w 



— [F(w , a^ + w d w F(w 0l ai 
w 



(4.31) 



W=WQ 



Since F is a rational function and F(wo, ai) is of order unity, the derivative d w F(wo, a{) 
cannot blow up. Therefore the second term in brackets can be ignored as wo — > 0, so using 
( |4.30| ) we arrive at the formula 



lim 

<li)->0 



da i 



ai 



-f(<ri)8 



((H) 


B\ Bi 







■f(<rt 



(4.32) 



which is valid under the assumptions on /, A and B given above. 

Now let us turn our attention to the soft gluon limit (|4.6| ). First we set the helicity 
of gluon 1 to +1 by setting r/i = in A n ^- This kills the % = 1 term in the third delta 
function in ( |1 . 1| ) . In the second delta function, the i = 1 term also vanishes trivially in the 
soft limit since A" — > 0. Particle number 1 therefore only appears in the integrals 



/ da ^ d ^ T- ( T7 ,8{\\-^A l )8{\\-i 1 B x ). 



(4.33) 
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The £1 integral is trivial and leads to 
1 1 



<r 2 (X\ 



1^2 



da i 



<Ji - a 2 (T\ - a n 



' W Ai-Ix 



(4.34) 



Now we are completely free to subtract from the argument of the delta function an amount 
which is equal to zero in the form of Xf/X} — Bi/Ai for any i ^ 2. (This is guaranteed to 
be zero by the Xf — £iP-* delta functions for i = 2.) Then we simply apply the formula 
Q4.32|) , once with i = 2 and once with % = n, to obtain the factor 



CT2 (X\ 



1^2 



a n -a 2 (X\y L (2 1> (nl) 



i \i 



AiA 2 AjA 



er, 



(n2> 



a 2 L(21>(nl) 



(4.35) 



The factor of 1/ (a n — 02) is needed to write the correct measure factor (|2.2|) for the (n — 1)- 
particle amplitude A(2, . . . , n). Gluon number one has now completely disappeared from 
the integral, leaving only the overall factor in brackets, in agreement with ( f4.6| ). 



4-4- The collinear limit 



First we consider the factor 



a 2 (X\ 



1^2 



da 1 



0"! 



, , Af Si 
M A["^ 



(4.36) 



which arises exactly as in the previous subsection. However, whereas we could there use 
( |4.32| ) for both i = 2 and i = n (since (1 2) and (1 n) were both going to zero), here we can 
only use ( |4.32j ) for i = 2 since only (12) — > in the collinear limit. Therefore the second 
term in brackets in ( |4.36j ) gives no contribution to the pole, and we only pick up the factor 



X\ 



&n ~ 02 A{ (2 1) a r 



0"2 



l-z 1 
* (21) 



(4.37) 



At this stage the integrals over the variables o\ and £1 associated with gluon number 
1 have been performed, but those associated with gluon 2 remain and we must rewrite 



the A2 dependence in terms of A = A2/VI — z - I n the ^2 integral this is accomplished by 
rescaling £2 in order to obtain 



d§2 
6 



8{X\-i 2 A 2 )8{Xl-i 2 B 2 ) 



1 



1 - z 
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/ ^^A 1 - &4 2 M(A 2 - &B 3 ). (4.38) 



The last delta functions to check are the ones of the form 

52 ( E ) = 52 + a? + E 



K i=l / \ i=3 



A 1 ~ ~ n ~ 

^A« + ^ 2 fc A« + E^ A i 

^ A 1 ~ ■ ~ ■ n 



(4.39) 



\-zA a 

i=3 

where in the first line we used the fact that we already integrated out £1 setting it to 
£1 = \\/Ai, in the second line we used the fact that we integrated out o\ setting o\ = 02, 
and in the third line we used the fact that A1A1 = ^ A 2 A2- Of course, we know that £2 
will eventually be set by a delta function to the value A2/A2, so we may as well write the 
final line as 

keeping in mind that A2 = yl — z\ and £2 = V^- ~ z ^2- 

What remains has precisely the structure of the amplitude A(p, 3, . . . , n), together 
with the extra factors in brackets from (|4.37|) and (f4.38|) , in complete agreement with the 



collinear limit (|4.7| ). The conjugate of this equation follows from the parity transformation 
discussed in section 4.2. The most notable fact following from that analysis is that the 
pole arises from the root satisfying a± — <J2 — [1 2] . One might attempt to prove the last 
collinear limit (|4.8|) by combining the above discussion with this observation. Then, the 
different z dependence might arise from the A dependence of the fermionic integrals. 



5. Conclusions and Speculations 

In this paper we have presented strong evidence that the formula ([lj]) encodes the 
complete tree-level ^-matrix of Yang-Mills theory in four dimensions. Explicit calculation 
has now shown that (|1.1|) agrees with YM theory for all MHV and MHV amplitudes, 
as well as all 6-particle non-MHV amplitudes. Moreover the analysis of section 4 shows 
that for any n, ( |TTT| ) satisfies a number of important properties required of gauge theory 
amplitudes, including parity symmetry. Many interesting directions remain open. 
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Of primary importance is to understand the connection between the formula ( |1 . 1|) , 
which was obtained in following the suggestion in |IJ] that one should consider a single 
instanton of degree d in the topological B-model on P 3 ! 4 , and the diagrammatic procedure 
of !|, in which arbitrary amplitudes are built out of d disconnected amplitudes, each of 
degree 1. We suspect that formulating a proof that ( |1.1|) factorizes correctly onto multi- 
particle poles would essentially amount to proving the equivalence of (|TTT|) and the rules 
of ||, simply because the factorization properties are completely manifest in the latter. 

The numerical coefficient in front of (|1 . 1| ) was fixed by comparing with gauge theory. 
We have not computed this coefficient independently in the B-model. It is conceivable 
that the degree d contribution and the separated degree 1 contributions (as well as other 
contributions) have to be added together to fully reproduce the normalization of the gauge 
theory scattering amplitudes. It is also possible that the B-model has some huge symmetry 
group which relates the connected instanton contribution ( p. . 1|) to the fully disconnected 
instantons of [[)]. 

Of course, even forgetting for the moment about the B-model, it would also be very 



interesting to prove rigorously that the formula (|1 . 1|) is the tree- level S'-matrix of Yang-Mills 
theory. To this end it would be useful to understand better the mathematical structure of 
the equations ( |2.6|) , and in particular to learn how many roots they have for general n and 
d (i.e., what is the degree of the corresponding Groebner basis). These numbers might be 
related to some interesting invariants of P 3 ' 4 or of its moduli space of holomorphic curves, 
and perhaps the equality of N n ^d and N n ^ n -d-2 could be understood in this language. 

Finally, all of our considerations have applied to the tree-level S-matrix in gauge 
theory. An obvious next step of great interest would be to see what light the topological 



B-model can shed on one-loop calculations [17 
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